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ABSTRACT 

We investigate the clustering properties of galaxies in the recently completed ELAIS- 
Sl redshift survey through their spatial two point autocorrelation function. We used 
a sub-sample of the ELAIS-Sl catalog covering approximately 4 deg^ and consisting 
of 148 objects selected at Ibpm with a flux > 0.5 mJy and redshift z < 0.5. We 
detected a positive signal in the correlation function that, in the range of separations 
1 — 10 /i“^Mpc is well approximated by a power law with a slope 7 = 1.4 ± 0.25 and 
a correlation length sq = 5.4 ± 1.2 /i“^Mpc, at the 90 % signihcance level. This result 
is in good agreement with the redshift-space correlation function measured in more 
local samples of mid infrared selected galaxies like the IRAS PSCz redshift survey. 
This suggests a lack of significant clustering evolution of infrared selected objects out 
to z = 0.5 that is further confirmed by the consistency found between the correlation 
functions measured in a local {z < 0.2) and a distant (0.2 < z < 0.5) subsample 
of ELAIS-Sl galaxies. We also confirm that optically selected galaxies in the local 
redshift surveys, especially those of the SDSS sample, are signihcantly more clustered 
than infrared objects. 
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1 INTRODUCTION 

Investigating the redshift-space distribution of galaxies has 
long been regarded as a fundamental aspect of observational 
cosmology. The primary statistical tool for characterizing 
galaxy clustering is the spatial two point correlation func¬ 
tion, ^(s) since, in the current paradigm of structure for¬ 
mation, the galaxy two-point correlation function is directly 
related to the initial power spectrum of mass density fluc¬ 
tuations. 

This is true as long as galaxies trace the underlying 
mass density held. However, it is now well established that 
the clustering of galaxies at low redshift depends on a vari¬ 
ety of factors, implying that not all types of galaxies can be 
regarded as unbiased mass tracers. The clustering of opti¬ 
cally selected galaxies has been found to depend on galaxy 
luminosity (Norberg et al. 2002 and reference therein), mor¬ 
phological and spectral type (Hermit et al. 1996, Zehavi et 
al. 2002, Magdwick et al. 2003). On the other hand, the 
clustering of infrared selected galaxies seems to depend on 


their infrared color (Hawkins et al. 2001) rather than on 
luminosity (Szapudi et al. 2000). 

The evidence that different galaxy populations might 
give different biased pictures of the mass distribution has 
complicated but also enriched the interpretation of galaxy 
clustering. Indeed, the very fact that the spatial clustering 
of galaxies is related to their physical properties represents 
an important observational test for all theories of galaxy for¬ 
mation. In particular, strong constraints on galaxy evolution 
models can be obtained by measuring the relative clustering 
of different extragalactic objects as a function of redshifts, 
for which very deep galaxy samples are required. 

Several deep redshift surveys of optically selected galax¬ 
ies, like the Keck surveys of the GOODS-north (Wirth et 
al. 2004) and GOODS-south (Vanzella et al. 2004) fields 
and the DEEP2 redshift survey (Davis et al. 2003), are cur¬ 
rently being performed. First results based on early data 
look very promising indeed. The analysis of Coil et al. (2004) 
has shown that the two-point correlation function of DEEP2 
galaxies with a median redshift 2 = 1.14 is consistent with 
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that measured by Adelbergeret al. (2003) in a very deep 
{z ~ 3) sample of Lyman break galaxies but is significantly 
smaller than the correlation measured in the local {z ~ 0) 
2dF galaxy sample. Since 2dF 6j-selected galaxies are known 
to trace the underlying mass density field in a unbiased way 
(Verde et al. 2002, Lahav et al. 2002), the smaller correla¬ 
tion measured in DEEP2 survey implies that this is not a 
strongly biased sample of objects either, since the clustering 
of the dark mass is expected to decrease with redshift in a 
similar way (Coil et al. 2004). 

Mid-infrared selected galaxies also constitutes a very in¬ 
teresting population of objects since they are also known to 
trace the underlying mass distribution in the local universe 
in an almost unbiased fashion. More precisely, it has been 
found that the fluctuations in their number density, 5g is 
related to the underlying mass overdensity field Sm through 
a simple, linear biasing relation 5g ~ 1.25m (Tegmark, Zal- 
darriaga & Hamilton 2001 and Taylor et al. 2002). Another 
reason why mid-infrared selection is interesting is that lu¬ 
minosity in this band is approximately proportional to star- 
formation rate (independent of dust), thus a mid-infrared 
sample of galaxies will highlight the distribution of star- 
formation activity at a particular epoch. Clearly, it would 
be very important to quantify the clustering evolution of 
mid-infrared selected objects and compare it with that of 
optically selected galaxies. This is indeed the main goal 
of this work in which we analyze the clustering properties 
of sample of mid-infrared selected galaxies extracted from 
the ELAIS redshift survey (Oliver et al. 2000) and extend¬ 
ing out to a redshift z = 0.5. According to La Franca et 
al. (2004), two main spectroscopic classes have been found 
to dominate the extragalactic population of these objects: 
star-forming galaxies (from absorbed to extreme starbursts: 

~ 10® — lO^^Lo), which account for 75% of the 
sources, and active galactic nuclei (excluded from this anal¬ 
ysis) which account for 25% of the sources. About 20% of 
the extragalactic ELAIS sources are dust-enshrouded star- 
burst galaxies, while passive galaxies are essentially absent 
from the sample. 

Our analysis is performed in redshift-space and thus 
complements the previous work of Gonzalez-Solares et al. . 
(2004) who measured the angular correlation properties of a 
similar sample of ELAIS galaxies that from which they have 
inferred their spatial correlation properties. 

The outline of this paper is as follows. In section 2 we 
describe the ELAIS-Sl galaxy redshift survey that we an¬ 
alyze in this work. In section 3 we discuss our method of 
estimating the two-point correlation function, assess its ro¬ 
bustness and evaluate its statistical uncertainties. The main 
results are presented in section 4, and discussed in section 
5, in which we also draw our main conclusions. 


2 THE ELAIS SI SAMPLE 

The European Large-Area ISO survey (ELAIS, Oliver et al. 
2000; Rowan-Robinson et al. 2004) is the largest Open Time 
programme conducted by the ISO satellite (Kessler et al. 
1996). It covers an area of 12 deg^, divided in four fields 
(Nl, N2 and N3 in the northern hemisphere and SI in the 
southern one ) distributed across the sky in order to de¬ 
crease the biases due to cosmic variance. The survey bands 


are at 6.7, 15, 90 and 170pm; the 15pm one presents the 
highest density of galaxies (Serjeant et al. 2000, Gruppioni 
et al. 2002, La Franca et al. 2004), making it the best choice 
for a study of their clustering properties. In this work we 
concentrate on the southern area, SI. 

This survey is made of nine raster observations, each 
covering 40 x 40 arcmin^. The final analysis catalog at 15 
pm in the SI field has been released by Lari et al. (2001) 
covering an area of 2 x 2 deg^ centered at a (2000) = 00^ 
34"* 44.4“, (5(2000) = -43° 28' 12". It includes 462 mid-IR 
sources down to a flux limit of 0.5 mjy. We have restricted 
our analysis to a highly reliable subsample of 406 objects 
(La Franca et al. 2004). 

A detailed description of the optical classification of the 
ELAIS-Sl sources, size and completeness function of the ar¬ 
eas used in our study, as well as the observed counts for 
each class of sources (normal galaxies, Starburst galaxies and 
AGN) are presented by La Franca et al. (2004). The mea¬ 
sure of the evolution of star-forming galaxies has been inves¬ 
tigated by Pozzi et al. (2004) and Gruppioni et al. (2005), 
while a first estimate of the luminosity function for type-1 
AGN has been presented by Matute et al. (2002) 

The central raster, S1_R5, has been observed three 
times and thus represents the deepest part of the sample. 
As a consequence the ELAIS-Sl area has been divided into 
two regions: a) the central and deepest part (S1_R5) which 
covers 0.55 deg^ and reaches a 20% completeness at fluxes of 
0.7 mJy, and b) the remaining area (Sl_Rest) which covers 
3.65 deg^ and reaches a 20% completeness at fluxes of 1.1 
mJy. 

The optical analysis of the 406 objects brings 332 (80%) 
optical identification on CGD exposures down to R~23, and 
290 spectroscopic classifications (90% of the optically iden¬ 
tified sample). Of these, 93 were found to be stars and 199 
extragalactic sources. Among these, we kept only the galax¬ 
ies with 2 < 0.5, because according to analyses based on 
the mean optical/mid-IR ratios of the sources, the sample 
is virtually 100% spectroscopically complete down to this 
redshift (La Franca et al. 2004). We have excluded all the 
AGN type sources. This left us with 148 redshift determined 
galaxies, 48 belonging to the deeper S1_R5 raster and the 
remaining 100 galaxies to the surrounding areas Sl_rest. 

The angular distribution of ELAIS-Sl galaxies in our 
sample is displayed in Fig.l (filled dots). The galaxy surface 
density is higher in the central region corresponding to the 
deeper S1_R5 raster. 

The continuous line histograms in both panels of Fig. 
2 shows the redshift distribution of the ELAIS-Sl galaxies 
in the S1_R5 (top) and SlJLest (bottom) rasters. Gonzalez- 
Solares et al. . (2004) have shown that the redshift distri¬ 
bution of ELAIS objects obtained from follow-up spectro¬ 
scopic observations and photometric redshifts is consistent 
with that predicted from the ELAIS luminosity function of 
Pozzi et al. (2004). This redshift distribution will used in 
the next section to estimate the galaxy correlation function 
and to construct the mock ELAIS catalogs. 


^ Data and related papers about the ELAIS southern survey are 
available at: http : //www.fis .uniromaS. it/ ~ ELAIS 
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Figure 1. The angular distribution of ELAIS-Sl sources. The 
solid square is the central raster S1_R5, observed three times. 


3 MEASURING THE REDSHIFT SPACE 
CORRELATION FUNCTION C(S') 

Our purpose is to estimate the redshift-space two point cor¬ 
relation function of ELAIS-Sl galaxies and its uncertainties. 
It is worth stressing here that clustering in redshift-space is 
systematically different from the true one in real-space. On 
small scale, velocity dispersion in virialized systems smears 
out structure along the line of sight, hence depressing the 
amplitude of ^(s) at small separation. On large scales, coher¬ 
ent infall motions in high density regions or outflows from 
low density environments enhance structure along the line 
of sight and thus increase the amplitude of f(s) on the cor¬ 
responding scales. Here we make no attempt of correcting 
for redshift-space distortions. Therefore, direct comparisons 
will only be made with similar measurements of the redshift- 
space two point correlation function of other galaxy redshift 
catalogs. 


3.1 Estimate of ^(s) 

Several different methods have been proposed to estimate 
the two point correlation function. In this work we use the 
Landy & Szalay (1993) [LS] estimator that, according to 
extensive numerical tests, turned out to be the best suited 
for astrophysical applications (Kerscher, Szapudi and Szalay 
2000 ). 

With the LS estimator a reliable and unbiased esti¬ 
mate of ^(s) is obtained by cross-correlating the real dataset 
with a sample of fake objects lacking any spatial correlation 
which, however, suffers from the same selection effects and 
incompleteness of the real catalog. To obtain a sample of 
synthetic objects we first create a random sample of points 
uniformly distributed over the same area as the real dataset, 
and then we generate their redshift distribution constrained 
to match that of the real objects determined by Pozzi et al. 



Figure 2. The redshift distribution of the sources in the S1_R5 
raster (top panel) and in Sl_Rest (bottom panel). The solid lines 
histogram indicates the redshift distribution of galaxies in the real 
ELAIS-Sl sample, while the dashed-line histogram indicates the 
average redshift distribution of objects in the 30 mock ELAIS-Sl 
catalogs. 


(2004). For each random catalog we generate 10 times the 
number of fake objects as real galaxies. 

The LS estimator is: 


, , ^ Ndd{s) _ r^ NpRis) 

’ Nrr{s) {npj Nrr{s) {npj^ ’ 


( 1 ) 


where Npp, Ndr and Nrr as the weighted data-data, data- 
random and random-random pair counts with redshift-space 
separation s. np and ur are the mean densities of galaxies 
and random objects, respectively. 

The flux limited nature of our dataset is accounted for 
by weighting each galaxy by 
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147vn{zi)J3{s)' 

where n{zi) is the observed space density of galaxies at red- 
shift Zi and Jsis) = ^{r)r^dr (Efstathiou 1988). To imple¬ 
ment this weighting scheme, that requires prior knowledge 
of ^(s), we have adopted an iterative procedure in which we 
first compute ^(s) using Wi = 1 and then we fit a power law 
to the results in the range 1 — 10 /i“^Mpc. Then we compute 
the new weights Wi by truncating the maximum value of J3 
to 1500. Only one iteration was needed to obtain a stable 
result. 

To check the robustness of our results we have varied 
both the estimator of ^(s) and the weighting scheme. More 
precisely, we have also used the estimators of Hewett (1982), 
Davis & Peebles (1983) and Hamilton (1993). Moreover, for 
the case of the LS estimator only, we have implemented two 
alternative weighting schemes: the case of no weight (i.e. 
Wi = 1) and the case in which we weight by the inverse of 
the selection function (i.e. Wi = nD/n{zi)). In all cases ex¬ 
plored we have found that the various estimates of ^(s) were 
consistent, within the errors, in the range 1 — 10 /i“^Mpc. 


3.2 Errors estimation 


To quantify the uncertainties in the estimate of ^(s) we have 
measured the correlation function in a sets of 30 independent 
mock galaxy catalogs designed to mimic the properties of the 
real ELAIS-Sl sample. 

These mock catalogs were constructed from the N-body 
numerical experiment labelled L3S performed by Cole et al. 
(1998). Their simulation assumes a flat ACDM model cos¬ 
mology with flm = 0.3, A = 0.7 an rms fluctuation of the 
mass contained in spheres of radius 8 /i“^Mpc, as = 1.13 and 
and a CDM power spectrum with shape parameter T = 0.25. 
The simulation box side is 345.6 /i“^Mpc and has 192® par¬ 
ticles. 

To construct each mock catalog we have used three dif¬ 
ferent output of the simulations (corresponding to redshifts 
z—0, z=0.33 and z=0.5) and took advantage of the peri¬ 
odic boundary conditions to obtain different configurations 
of points in each box by re-centering the coordinate system 
on a particle chosen at random. We have then stacked six 
boxes in order of increasing redshift (three of them at z=0, 
two at z=0.33 and one at z=0.5), placed the observer at the 
center of the first box and identified the X coordinate with 
the direction of stacking. We have then identified the direc¬ 
tion of the X axis with the center of the ELAIS field and 
considered all particles within an area of 2 x 2 deg^ from the 
field centre. To assign each particle a redshift we have used 
the distance to redshift relation for a flat, ACDM universe 


r = 



dz! 

zt)^ -\- Ha 


( 3 ) 


where r is the comoving distance, and then added the line 
of sight component of the particle peculiar velocity. Finally, 
a population of mock galaxies have been extracted from the 
particles through a Montecarlo rejection procedure designed 
to match the observed redshift distribution of real ELAIS 
galaxies in both the inner 40 x 40 arcmin^ Sl_R5-like, and 
the outer Sl_Rest-like samples. This procedure has been re¬ 
peated to obtain 30 mock ELAIS-Sl samples, containing 



Figure 3. Average ^(s) measured in the 30 mock ELAIS-Sl cata¬ 
logs (filled dots) and the rms scatter around the mean (errorbars). 
Open dots shows the “true” ^(s) of all particles in the simulation 
contained in a region with the same volume of the ELAIS-Sl 
sample. The straight line represents the best power law fit to the 
“true” ^{s) in the range 1 — 10 /i“^Mpc. The values of the best 
fit parameters are also shown. 


^ 120 objects each, for which we have evaluated the two- 
point correlation function. The mean redshift distribution of 
fake objects in the 30 catalogs is shown in of Fig. 2 (dashed 
line histogram) for the S1_R5 (top panel) and Sl_Rest (bot¬ 
tom) rasters 

The filled dots in Fig. 3 shows the average ^(s) in the 
30 mock catalogs. The errorbars represent the rms scatter 
around the mean values. They constitutes our estimate of er¬ 
rors in the measurement of two point correlation function of 
the real ELAIS-Sl galaxies These errors account for cosmic 
variance and sample noise, that constitutes the main source 
of uncertainty. The large size of the sample guarantees that 
the ’integral constraint’ correction is only 5% (Gonzalez- 
Solares et al. 2004) and thus can be neglected in the error 
budget. 

The open dots in Fig. 3 shows the “true” ^(s) of the 
simulation, measured by considering all particles in one of 
our ELAIS-like samples. The “true” ^(s) is consistent with 
the average ^(s) of the mocks, indicating that our method 
for estimating ^(s) is indeed unbiased. 

The mock ELAIS-Sl catalogs mimic the geometry and 
selection effects of the real sample and account for the spa¬ 
tial clustering and its evolution in a flat ACDM universe. 
However, they are not guaranteed to reproduce the cluster¬ 
ing properties of the ISO galaxies unless, of course, these 
trace the underlying mass density field out to 2; = 0.5. 

As we have verified a posteriori this is indeed the case, 
in the sense that the “true” ^(s) in the range 1 — 10 /i~®Mpc 
is well approximated by a power law whose best fit param¬ 
eters, displayed in Fig. 3 are consistent within the errors 
with those determined in the analysis of the real ELAIS-Sl 
sample presented in the following section. 
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Figure 4. The redshift-space correlation function for the ELAIS- 
S1 galaxy sample (filled dots). Error bars have been calculated 
using the 30 ELAIS-like mock catalogs. The solid line represents 
the power law best fit ^(s) = (s/sq) > to the data in the range 
1 — 10 /i“^Mpc. The best fit parameters are indicated in the plot. 

4 RESULTS 

The filled dots plotted in Fig. 4 represent the two-point cor¬ 
relation function of the ELAIS-Sl galaxies computed us¬ 
ing the LS estimator. The errorbars, evaluated from the 
30 ELAIS-Sl mock catalogs, are the same shown in Eig. 
3. Clearly, ^(s) is well approximated by a power law out to 
separations of 10 /i“^Mpc. In the range 1 — 10 /i“^Mpc the 
best fit power law model ^(s) = {s/so)~^, has a correlation 
length of So = 5.4± 1.2 /i“^Mpc and a slope 7 = 1.45±0.25, 
both determined at the 90% confidence level. Including the 
correlation of galaxy pairs at separation s = 0.6 /i“^Mpc, 
also shown in the figure, does not modify this result appre¬ 
ciably. 

Breaking down the sample by redshift does not change 
results significantly either. Indeed, in the local sample com¬ 
posed by 82 objects at z < 0.2 ^(s) is still well approxi¬ 
mated by a power law in the range 1 — 10 /i“^Mpc with best 
fit parameters so = 5.4 ± 1.6 /i“^Mpc and 7 = 1.6 ± 0.4. 
These values agree with those found for the distant sample 
of 66 objects at at 0.2 ^ 2 ; < 0.5 for which we have found 
So = 5.1 ± 1.6 /i“^Mpc and 7 = 1.4 ± 0.4, 

Finally, we have verified that evaluating errors from 100 
bootstrap realizations of the ELAIS-Sl sample rather than 
from the mock catalogss does not change significantly the 
results as the correlation length only decreases by ^ 10 % 
and the slope becomes ~ 4% flatter. 


5 DISCUSSION AND CONCLUSIONS 

In this work we have evaluated the redshift-space two point 
correlation function of mid-infrared selected galaxies in the 
deep (2 0.5) ELAIS-Sl catalog. We found a significant, 

positive correlation signal at separations ^ 10 /i“^Mpc, 


where the two point correlation function is well approx¬ 
imated by a power law model ^(s) = {s/so)~^ with a 
correlation length of so = 5.4±1.2 /i“^Mpc and a slope 
7 = 1.45 ± 0.25, with errorbars referring to a 90% confi¬ 
dence level. These results have been obtained using the LS 
estimator for ^(s) and by evaluating the errors from a set of 
30 mock ELAIS-Sl catalogs. These results are robust, in the 
sense that they do not change significantly when varying the 
method of estimating ^(s) (see section 3.1) or when using 
different strategies to assess the errors (sections 3.2 and 4). 

It is interesting to compare our results with those 
obtained from similar analysis of galaxy clustering in 
redshift-space. Table 1 shows our result together with 
the corresponding results obtained from some of the ma¬ 
jor galaxy redshift surveys, characterized by their pass- 
band/wavelength of selection (column 2 ) and the redshift 
ranges they cover (column 3). It is worth stressing that de¬ 
viations of the two point correlation function from a pure 
power law shape are more serious in redshift-space than in 
real-space, making it difficult to compare results obtained 
from the analyses of different galaxy samples. To ensure a 
fair comparison, all best fit parameters listed in columns 4 
and 5 of Table 1 refer to ranges of separations where all 
measured ^(s) are well approximated by a power law. These 
ranges turned out to be very close to that of [1 —10] /i“^Mpc 
considered in our analysis, although some of the parameters 
in the Table have been obtained by pushing the estimate 
of ^(s) down to scale as small as 0.1 /i“^Mpc (the 2dFGRS 
case) or up to separations as large as 16.4 /i“^Mpc (as in 
the case of the LCRS sample). 

Our results are fully consistent with those obtained from 
the analysis of the PSC 2 survey (Hawkins et al. 2001), that 
consists of ~ 15, 000 IRAS galaxies selected at 60/im, i.e. in 
a mid-infrared band similar to that of ELAIS galaxies. This 
sample is, however, much more local than ours and thus 
the agreement between the two results indicates that the 
clustering properties of mid-infrared selected objects do not 
evolve significantly between 2 = 0 and 2 = 0.3. This conclu¬ 
sion is corroborated by the consistency found between the 
two measurements of ^(s) performed in the local (2 < 0 . 2 ) 
and distant (0.2 ^ 2 < 0.5) ELAIS-Sl subsamples. Unfortu¬ 
nately, the large uncertainties in our estimate of ^(s), that 
mainly result from the sparseness of the ELAIS-Sl galaxy 
catalog, do not allow to set strong constraints on the clus¬ 
tering evolution. 

The lack of significant evolution in the ELAIS galaxy 
clustering has already been noticed by Gonzalez-Solares et 
al. (2004). In their analysis, that consisted in deprojecting 
the angular correlation function of ELAIS galaxies via Lim¬ 
ber equation, they have measured a real-space two point cor¬ 
relation function with a slope (7 = 2.04 ±0.18) and a corre¬ 
lation length (ro = 4 . 3 ('lg '7 /i“^Mpc) that are fully consistent 
with those measured in the PSC 2 catalog (7 = 2.04 ± 0.18 
and ro 3.7 /i“^Mpc, Jing, Borner & Suto (2002)). It is 
worth stressing that difference between our result and that 
of Gonzalez-Solares et al. (2004) originates from systematic 
redshift space distortions that affect our analysis and result 
in a shallower slope and a larger correlation length of the 
two point-correlation function. 

hocusing on the mid-infrared objects is of considerable 
interest since they trace the underlying mass distribution 
in the local universe. The consistency that we have found 
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Table 1. Clustering of Different Galaxy Redshift Surveys: the Parameters of the Power Law Model 


Survey 

Selection Band 

Az 

so( b '^Mpc) 

7 

ELAIS-Sfi 

15/im 

0.0-0.5 

5.40 ± 1.20 

1.45 ±0.25 

PSC^2 

60/im 

0.004-0.1 

4.77 ±0.20 

1.30 ±0.04 

CfA23 

B-band 

0.0-0.05 

~ 7.5 

~ 1.6 

ORS'^ 

B-band 

0.0-0.027 

7.60 ± 1.20 

1.60 ±0.10 

LGRS® 

R-band 

0.033-0.15 

6.3 ±0.3 

1.86 ±0.03 

SDSSS 

r*-band 

0.019-0.13 

8.0 

~ 1.2 

2dFGRS® 

bj-band 

0.01-0.20 

6.82 ±0.28 

1.57 ±0.07 


^This analysis 
^Hawkins et al. (2001) 

®de Lapparent, Geller &; Huchra (1988) 
“^Hermit et al. (1996) 

^Tucker et al. (1997) 

®Zehavi et al. (2002) 

^Hawkins et al. (2003) 


between the ^(s) in the real sample and the one measured 
in our mock ELAIS-Sl catalogs in which galaxies are iden¬ 
tified with the particles of the parent N-body simulation, 
suggests that mid-infrared selected galaxies are indeed al¬ 
most unbiased tracers of the underlying mass density field 
out to 2 = 0.5. 

It is interesting to compare our results with those 
obtained using optically-selected galaxy samples, some of 
which is listed in Table 1. We confirm the well known fact 
that infrared selected objects are significantly less clustered 
than optical galaxies, although the discrepancy depends on 
the depth of the sample and (mainly) on the selection band, 
being larger for SDSS, r* selected objects than for 2dF 
galaxies selected in the bj band. 

Investigating how this discrepancy change with redshift, 
i.e. how the relative clustering of different objects evolve 
with time, is of fundamental importance since it will allow to 
cast light on the galaxy formation and evolution processes. 
While this work constitutes an early step in this direction, 
strong observational constraints for theories of galaxy evolu¬ 
tion will only be obtained by extending and improving sam¬ 
ples like ELAIS-Sl, which implies carrying out larger and 
deeper infrared redshift surveys like that already planned in 
the SWIRE programme (Lonsdale et al. 2003). 
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